Center for Turbulence Research 
Annual Research Briefs 1992 



Plane mixing layer vortical structure kinematics 

By R. L. LeBoeuf 


The objective of the current project was to experimentally investigate the struc- 
ture and dynamics of the streamwise vprticity in a plane mixing layer. The first 
part of this research program was intended to clarify whether the observed decrease 
in mean streamwise vorticity in the far-field of mixing layers (Bell & Mehta 1992) 
is due primarily to the “smearing” caused by vortex meander or to diffusion. Two- 
point velocity correlation measurements have been used to show that there is little 
spanwise meander of the large-scale streamwise vortical structure. The correlation 
measurements also indicate a large degree of transverse meander of the streamwise 
vorticity which is not surprising since the streamwise vorticity exists in the inclined 
braid region between the spanwise vortex core regions. The streamwise convection 
of the braid region thereby introduces an apparent transverse meander into mea- 
surements using stationary probes. These results were corroborated with estimated 
secondary velocity profiles in which the streamwise vorticity produces a signature 
which was tracked in time. 

1. Motivation and objectives 

An extensive data set consisting of single-point mean and turbulence statistics 
has been obtained for a two-stream mixing layer (Bell & Mehta 1989b, 1992). The 
plane unforced mixing layer originating from laminar boundary layers was examined 
in order to quantify the development of streamwise vorticity which previously was 
identified only through flow visualization studies (e.g. Bernal & Roshko, 1986). The 
mean streamwise vorticity derived from the mean velocity field shows a continuous 
decrease in magnitude with streamwise distance from its nearfield occurrence. It is 
unclear whether the decrease in mean vorticity is a result of diffusion of the stream- 
wise vorticity or due to meander of concentrated vorticity. Based on comparisons 
with forced streamwise vortex meander in a boundary layer, Bell & Mehta (1992) 
argued that the observed decrease of the mean vorticity in the far-field mixing layer 
was more likely a result of diffusion. 

Townsend (1976) showed that the governing equations for a free- shear flow admit 
to self-preserving solutions for sufficiently high Reynolds numbers. The resulting 
“self-similar” mean and Reynolds stress profiles become functions of single length 
and velocity scales. Previous measurements (Bell & Mehta 1992) have indicated 
that the streamwise vorticity persists even in what would normally be considered 
the “self-similar” region (where a linear mixing layer growth rate and asymptotic 
peak Reynolds stresses were achieved). The peak streamwise vorticity and the 
secondary shear stress (uu>), which was strongly correlated with the streamwise 
vorticity, were found to exhibit significant levels in this region (although they de- 
creased with streamwise distance to levels comparable with the noise threshold). It 
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is important for the establishment of the criteria for “self- similarity” to investigate 
whether the measured decay is due to true diffusion of the streamwise vorticity or 
an artifact of meander. In addition, this assessment will have important implica- 
tions regarding the ability of the layer to enhance mixing and reaction rates in the 
far-field. 

To resolve the questions regarding the persistence of streamwise vorticity in the 
far-held, it was proposed to perform two-point cross-wire measurements of the ve- 
locity held (Bell 1990). The dependence of the velocity cross-correlation on the 
hxed probe location is considered a good indicator of the stationarity of the stream- 
wise vortex location. Additional information regarding meander can be obtained 
from instantaneous velocity prohles. These were estimated using a newly developed 
technique and the current apparatus (LeBoeuf & Mehta 1992). The facility and 
how conditions that were used for the current study are similar to those of Bell & 
Mehta (1992). 

2. Accomplishments 

2.1 Experiment and apparatus 

The experiments were conducted in a Mixing Layer Wind Tunnel specifically 
designed for free-shear how experiments (Bell & Mehta 1989b). The wind tunnel 
consists of two separate legs which are driven independently by centrifugal blowers 
connected to variable speed motors. The two streams merge at the sharp trailing 
edge of a slowly tapering splitter plate; the included angle at the splitter plate edge, 
which extends 15 cm into the test section, is about 1°. The test section is 36 cm in 
the cross-stream direction, 91 cm in the spanwise direction, and 366 cm in length. 
One side-wall is adjustable for streamwise pressure gradient control and slotted for 
probe access. 

In the present experiments, the two sides of the mixing layer were set to 15 m/s 
and 9 m/s for a velocity ratio, r = U 2 /U 1 = 0.6 [A = (J7j — U 2 )/(U\ + U 2 ) = 
0.25]. For these operating conditions, the measured streamwise turbulence level 
{u 9 JUt) was approximately 0.15% and the transverse levels ( 1 / jU € and w , /U e ) were 
approximately 0.05%. The mean core-how was found to be uniform to within 0.5% 
and cross-how angles were less than 0.25° (Bell &; Mehta 1989a). The boundary 
layers on the splitter plate were laminar at these operating conditions. 

Measurements were made using two independently traversed cross-wire probes. 
The probes could be rotated in order to measure how in two-coordinate planes. 
The geometry of the instrumentation resulted in a minimum probe spacing of 7 
mm. One probe was mounted on a 2-D traverse which was designed and con- 
structed for the current work. The new traverse was manually controlled using 
an indexing stepper motor controller. The second probe was mounted on a pre- 
existing computer-controlled 3-D traverse. Both cross- wires were linked to a fully 
automated data acquisition and reduction system controlled by a DEC MicroVax II 
computer. The software required for multiple-probe measurements was developed 
for the current study. The Dantec cross- wire probes (Model 55P51) consisted of 
5 fim platinum-plated tungsten sensing elements approximately 1 mm long with 
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FIGURE 1. Mean streamwise vorticity (fix) at X = 57 cm, starting from -0.08, 
incremented 0.02. 


approximately 1 mm spacing. The probes were calibrated statically in the potential 
core of the flow assuming a ‘cosine-law’ response to yaw, with the effective angle de- 
termined by calibration. The calibration locations were selected to insure consistent 
angularity of the measurements of the two probes. The analog signals generated by 
NASA-built anemometers were DC offset and amplified (x 8) before being fed into a 
computer interface. The interface contained a fast sample-and-hold A/D converter 
with 16 bit resolution and a multiplexer for connection to the computer (Bell &c 
Mehta 1992). Individual statistics were averaged over 8192 samples obtained at a 
rate of 1500 samples per second. 

Data were obtained in several Y-Z planes between X = 17 and X = 128 cm 
with two orientations of the probe {uv and u w ) in order to identify mean vortex 
locations. This method yielded all three components of mean velocity, five indepen- 
dent components of the Reynolds stress tensor, and selected higher order products. 
Typically, 700-1200 points were measured on a rectangular cross-plane grid at a 
given streamwise station. The grid spacing (same in both directions, Y and Z ) was 
varied from 0.25-0.5 cm for upstream to downstream stations. The mean spanwise 
velocity (W) measurements were corrected for mean streamwise velocity gradient 
( dU/dY ) effects (see Bell and Mehta 1989b for details). The streamwise component 
of mean vorticity was computed using a central difference numerical differention of 
the measurements of the velocity components V and W. 

An example of the resulting measurement of the mean streamwise vorticity (fix) 
at X = 57 cm downstream of the splitter plate trailing edge is given in Fig. 1. The 
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selection of the traverse grid for the correlation measurements is also noted on the 
figure. Except for the decreasing peak magnitude and increasing size (for increasing 
streamwise distance), the nominally linear array of counter-rotating vortices was 
typical of all measurement locations except the first. At the first measurement 
location, X = 17, the streamwise vorticity had not yet organized into a single 
spanwise row. The correlation locations were nominally selected to cross the location 
of maximum mean streamwise vorticity. 

The two-point cross-correlation measurements required for the estimation of ve- 
locity profiles were acquired on grids of equally spaced probe positions. The spacing 
of measurement locations in the grids was 0.4 cm. The minimum distance between 
estimation and reference locations was limited by probe geometry to two grid spac- 
ings (i.e. 0.8 cm). Therefore, the estimated velocity profile is twice as dense as it 
could be using probes of the same geometry to measure at all grid locations simulta- 
neously. The grids were positioned 78 cm downstream of the splitter plate trailing 
edge for the measurements discussed in this paper. It is estimated that this station 
is located after the third pairing of the spanwise vortex structures. 

2.2 Results and discussion 

The transverse ( v ) and spanwise (w) velocity correlations in the spanwise ( Z ) 
and transverse (Y) directions, respectively, are shown in Figs. 2a and b for the 
streamwise location X = 57 cm from the splitter plate trailing edge. One noteworthy 
feature of the sets of curves for each station is the relative phase of each curve. 
There is generally a shift of the curves from right to left with increasing stationary 
probe location. There would be a one-to-one correspondence between the change 
in the probe spacing at the correlation extremes and the change in position of 
the stationary probe if no meander of the large scale structure is occurring. The 
transverse velocity correlations on a spanwise measurement grid approximate this 
no-meander relative phase shift. For example, Figs. 3a and 3b show the probe 
separation at the first minimum and first maximum (excluding zero spacing) of all 
of the correlation measurements versus stationary probe location. The slope of the 
resulting curves for the transverse velocity measurements (Fig. 3a) is very close 
to -1, indicating that meander in the spanwise (z) direction is negligible at this 
streamwise location. 

The spanwise velocity correlations across the mixing layer have a less consistent 
phase shift for changing stationary probe location. Typically, the shift of the first 
minimum is a fraction of the shift of the stationary probe for two sets of correlations 
(see Fig. 2b). This is reflected in Fig. 3b where the slope of the curve of first 
minimum versus stationary probe location for the spanwise velocity measurements 
has a much more gradual slope, which would indicate a large degree of apparent 
transverse meander. The slope of the curves may be linked to the aspect ratio of 
the mean streamwise vortices. 

A simple, more direct way to identify meander of the large-scale structure is to 
track them in temporal multi point velocity records. The presence of streamwise 
structure should appear clearly in profiles of the secondary velocity. In particu- 
lar, spanwise meander should produce a jittered signature in spanwise profiles of 
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FIGURE 2. Two-point cross- correlation at X = 57 cm. (a) Transverse velocity (u) 

in the spanwise ( Z ) direction at Y = 0 and Z = -8.5 cm ; -8.15 cm ; 

-7.8 cm ; -7.45 cm , (b) Spanwise velocity (u>) in the transverse (F) 

direction at Z = -5 cm and Y = -2.1 cm ; -1.4 cm ; -0.7 cm , 0.0 


cm 



FIGURE 3. Probe separation at extremes of correlations versus relative stationary 

probe location for X — 27 cm ; 47 cm ; 57 cm ; 78 cm ; 88 

cm ; 128 cm ; Line with slope equal to -1 . (a) Transverse velocity 

(v) in the spanwise ( Z ) direction, (b) Spanwise velocity (u>) in the transverse (y) 
direction. 
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transverse velocity and transverse meander should produce a jittered signature in 
transverse profiles of span wise velocity. Profiles of simultaneous velocity records 
with sufficient spatial resolution would be very difficult using cross- wire probes and 
would entail a significant increase in hardware costs. Use of optical techniques 
such as digital particle image velocimetry would probably be the best suited to this 
problem, but would also be costly. An alternative approach whereby the velocity 
profiles were estimated using one or two reference velocities was used to produce 
simultaneous velocity profiles. This approach required no further investment in 
hardware. 

The linear estimate of the velocity component u,(x') based on m reference (mea- 
sured) velocity components u ; (x (l) ) at n locations at the same time can be expressed 

as: „ _ 

n in 

««oo = ^2 ^2 Aij k uj(x ik) ). (i) 

*=1 j=l 

Minimizing the mean-square error, [t£^(x') — Uj(x')] 2 , of the estimate yields an equa- 
tion for the linear estimation coefficients 


n m 

EE Aij k Uj(xW)ui(xM) = it/(x(P>)ui(x'). (2) 

fc=l ;= 1 

Thus, given the cross-correlation tensors, the linear mean-square (MS) estimation 
coefficients can be determined. It can be shown [using Eq. (1)] that the left hand 
side of Eq. (2) is equal to ui(x(p))uj(x'). Hence, the estimation-reference cross- 
correlation tensor is equal to the measurements used to produce the coefficients. 

The mean-squared estimation error increases with the distance away from the 
measurement location (approaching u) ) , as the velocity loses correlation with the 
reference (Adrian & Moin 1988). The estimated signal can be amplified in a ratio- 
nal manner by including the velocity covariance at the estimation locations in the 
scheme to calculate the estimation coefficients. This optimization which includes the 
estimation covariance in addition to the two-point cross-correlation tensors forms 
the basis for the proposed technique. 

Consider the error vector, 

F — ( (3) 

V Uj(x') tti(x') - uj(x')ui(x') ) 

Then the estimation coefficients ( A i]k ) can be op timally determ ined with respect to 
both the estimation-reference cross-correlation, Uj(x')ui(x( k '>), and the estimation 
covariance, Uj(x')u,(x'), by minimizing ( F T F ), where the superscript T indicates 
transpose. The statistics of the estimates can be calculated for any value of the 
coefficients using known reference signals and estimation signals generated using 
Eq. (1). This optimization can be accelerated, however, by calculat ing the statistics 
of the estimates in terms of the reference cross-correlation tensor, u ; (x(*>)uj(x(p)), 
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and the estimation-reference cross- correlation tensor, u/(x(jd)u<(x'). By operating 
on Eq. (1), it can easily be shown that the estimation- reference correlation can be 
expressed as: 

71 771 

«j(*')“i(* ( *>) = 5Z 5Z i4 i»“«(* (,,) ) tt i(* ( * ) ) ( 4 ) 

p=l <7=1 

and the estimated Reynolds stresses can be expressed as: 

n 771 71 771 

u j (x')u,(x l ) = 53 /L S Aj„Aii k u q (xW)u,(xi k >). (5) 

p= 1 <7=1 k= 1 f=l 


Thus, the new technique requires the same cross- correlation tensor data a s the 
traditional mean-square estimate along with the estimation covariance data. The 
optimization described above is a non-linear, least-squares optimization which can 
be solved using standard techniques such as, for example, the modified Levenberg- 
Marquardt algorithm implemented in the IMSL subroutine DUNLSF. 

Another alternative to minimizing the mean-squared error is to maximize the 
correlation coefficient of the estimated signal and the true signal at the estimation 
location. When this criterion is used to determine the linear estimation coefficients 
for a two-component one-point estimate, the maximum correlation coefficient is 
independent of the energy of the estimate: 


Pit = 


( -uf(x' 

)a 



r 

2 
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where 








2u(xQu2(x^)ui(xW)l£2(x( fe )) 

fi(x')ui(xW) 
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and 


6 = 


2u 2 (x i ) 


ri(x')ui(xW) 




u(x')l/2(xW)Uj(xW) 

ti(x')ui(xW) 


( 3 ) 


Hence, once the maximum correlation coefficient has been calculated, the estimation 
coefficients can then be determined from 


A 2 — 


-Pub 
2 a 


(9) 


and 

A _ /),;u 2 (x') - A 2 u(x , )M 2 (^ (t) ) 


( 10 ) 
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consecutively. This restricts the energy of the estimate to the measured energy at 
the estimate location. 

Verification of the MS estimation and the new estimation schemes outlined above 
are detailed below. The verifications include comparisons of partial kinetic energy 
and the estimation-measurement correlation coefficient. Due to the simplicity of 
the linear estimation, these quantities (and many others) can be calculated using 
only the statistics required for the coefficient calculation (i.e. the time series do 
not have to be retained). For example, the partial kinetic energy of the estimated 
velocity signals can be expressed in terms of the Reynolds stresses given by Eq. (5). 
Similarly, the estimation-measurement correlation coefficient 

Uj(x')ui(x') /jj) 

P * [fi^(x l )u?(* , )J 1 / 2 

can be expressed in terms of the optimization data, u?(x') and u g (x^)u t (x^ k ^) y 
since the estimation-measurement correlation is given by: 

n m 

Uj(x')Ui(x') = ^2^2 Aj qp u 9 (xM)ui(xW) (12) 

p=l 9=1 


and uj(x') can be obtained from Eq. (5). By not requiring the calculation of such 
statistics from entire time series, optimum reference locations can quickly be de- 
termined without extensive calculation or data storage costs. Furthermore, the 
appropriate number of references to be used for time-series estimation can be se- 
lected using only two-point measurements. 

Cole et ai (1992) suggested comparison of the partial kinetic energy (uj + u]) in 
order to quantify the effectiveness of the estimation. The results of such a compar- 
ison for two sets of reference locations are shown in Figs. 4a and b. Note that in 
all figures, the results of the estimation schemes have been designated MS for the 
traditional mean-square estimate, LS for the new technique requiring the non-linear 
least-squares fit, and OP for the new technique which maximizes the estimation- 
measurement correlation coefficient. By taking two reference locations near the 
center of the mixing layer (Fig. 4a), the traditional stochastic estimate appears to 
recover a large fraction of the measured energy. This result is consistent with the 
findings for the jet mixing layer (Cole et a l. 1992). Unfortunately, most of the 
energy contribution is from the references and little energy is represented by the es- 
timates. It may also be noted in Fig. 4a that including the estimation covariance to 
the optimization for this choice of reference locations improves the energy recovery 
of the estimate very little (~ 10%). The result of taking reference locations near the 
edges of the mixing layer is shown in Fig. 4b. The traditional stochastic estimate 
for these choices of reference locations yield a negligible representation of the energy 
away from those locations. The estimation proposed here gave a very good repre- 
sentation of the energy; even better than that obtained for more centrally located 
references (as in Fig. 4a). The reason for the large improvement over mean-square 
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FIGURE 4. Transverse partial kinetic energy comparisons, (a) Two-point estimates 
with references at Y = -0.6 cm and Y = 0.2 cm and one-point estimates with 
reference at Y = -0.6 cm. (b) Two-point estimates with references at Y = -5.0 cm 
and Y = 3.4 cm and one-point estimates with reference at Y = 3.4 cm, Measured 

o ; Measured reference • ; Two-point MS ; Two point LS ; One-point 

MS — ; One-point LS ; One-point OP . 

estimation for this choice of reference locations is that the estimate-reference cross- 
correlation is much less than the estimation covariance. Therefore, minimizing the 
error in both quantities results in a better fit to the estimation covariance (which 
includes Uj and u This trade-off does not occur for reference locations near the 
center of the mixing layer. Consequently, the energy recovery for spanwise profiles 
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FIGURE 5. Transverse estimation-measurement correlation coefficients. Two-point 
estimates with references at Y = -5.0 cm and Y = 3.4 cm and one-point estimates 
with reference at Y = 3.4 cm. (a) Streamwise component ( u ). (b) Transverse 

component (v), Reference • ; Two-point MS ; Two-point LS ; One- 

point MS ; One-point LS ; One-point OP . 

near the center of the mixing layer was very poor for both estimation techniques. 
By taking spanwise grids near the mixing layer edge, a very good representation of 
the energy by the new estimation was obtained. 

Another measure of the success of an estimation is the correlation coefficient 
for corresponding estimated and measured velocities at the same location. If an 
estimated velocity matches the measurement exactly, then the correlation coefficient 
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FIGURE 6. Transverse velocity (v) measurement vs. estimate time series compar- 
isons for a reference at Y = 3.4 cm. (a) Y = 2.6 cm. (b)Y ~ -0.6 cm. (c) Y = 
-5.0 cm, Measured ; MS ; LS ; OP . 
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FIGURE 7. Estimated velocity profiles, (a) Transverse secondary velocity vectors 
at Z = -7.0 cm for a reference at Y = -3.4 cm. (b) Spanwise profile of transverse 
velocity at Y = 3.4 cm for references at Z = -5.0 cm and Z = -7.0 cm, starting 
from -0.1 m/s, incremented 0.08 m/s. 

would obviously be equal to one. The results of a calculation of the correlation 
coefficient for a high speed edge reference are shown in Figs. 5a and b. A high 
correlation between the estimation and corresponding measurement exists for most 
of the high-speed side locations. This is reflected in the corresponding time trace 
comparisons of Fig. 6a. It is no surprise that the low frequency trigger provided 
by the flow near the mixing layer edge does not yield estimates which are highly 
correlated with the high frequency signals near the mixing layer center (see Fig. 6b). 
The low frequency estimate resulting from the technique proposed here is largely an 
amplified version of the MS estimate. The amplification is a result of an attempt (by 
the optimization) to match the covariance of the estimate to that of the measured 
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signal. The estimation of flow on the opposite side of the layer is remarkably good 
as suggested by the high correlation coefficient for those locations (Figs 5a and b). 
Presumably, the large-scale structure providing the trigger at the high-speed edge 
produces similar signatures at the low speed edge, resulting in a good estimation of 
the flow there (Fig. 6 c). The streamwise component time series comparisons were 
left out for brevity since they were qualitatively similar to the transverse component 
comparisons of Figs. 6 a through 6 c. 

The results of estimating the one-component simultaneous multi-point velocity 
profile using the OP method is presented as Fig. 7a. This profile was generated for 
a transverse grid at Z = -7 cm using a reference at Y = -5 cm. Taylor’s hypoth- 
esis was used to approximate the spatial development from the temporal signals. 
This entailed calculating the relative spatial locations of velocity vectors as U c / fs, 
where f s is the sampling frequency and U c is the mean convection velocity. The 
secondary velocity vector plot clearly shows the spanwise rollers passing through 
the measurement station. In spite of the poor representation of the broad-band 
signals near the center of the mixing layer, the large-scale structure is adequately 
represented when edge references are selected. In particular, the relationship be- 
tween the kinematics of the large-scale streamwise and spanwise flow structures can 
be inferred from simultaneous estimates of velocity profiles on spanwise and trans- 
verse grids. For example, using two cross- wire probes makes it possible to estimate 
three velocity components in one profile (in one direction) and two velocity compo- 
nents in another profile (in another direction) provided the grids have a common 
reference location. This can be achieved by measuring reference signals U\ and u<i 
(or us) at the common grid location and simultaneously measuring U\ u 3 (or U2) 
in the grid for which three velocity components are desired. With the use of such 
a combination of grid locations, it was possible to estimate the transverse velocity 
(t>) as shown in Fig. 7b. The signature of the streamwise vortices which ride the 
braid between spanwise vortices appear each time the spanwise vortices carry them 
through the measurement grid. This gives rise to the juxtaposed spanwise sets of 
islands in Fig. 7b between time steps -17 and -12. The lack of significant structure 
to structure spanwise relative motion of this signature observed in longer simulta- 
neous time series profiles is indicative of the low degree of spanwise meander which 
was also inferred from the correlation measurements. 

3. Future plans 

During the remainder of the program, phase locked velocity measurements will 
be obtained by acoustically forcing the mixing layer. The hardware required for 
this phase of the program has already been prepared. By establishing a periodic 
sequence of spanwise vortices, it should be possible to examine the detailed structure 
of the streamwise vortices. By measuring the scale and spacing of the streamwise 
vorticity in the braid region of the mixing layer vortices, more direct comparisons of 
the wind tunnel measurements can be made with the direct numerical simulations 
of the same flow (Rogers, M. M. &; Moser, R. D. private communication). In 
particular, initial conditions (streamwise vortex strength and distribution) will be 
measured in detail so that they may be used as input to numerical simulations. 
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During the latter part of the program, some scalar mixing studies using heat to tag 

one of the streams will also be conducted. An extension of the estimation to include 

space-time correlations will also be attempted if time permits. 
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